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■ Abstract 

We show herein that a pattern based on FGLM techniques can be used 
for computing Grobner bases, or related structures, associated to linear 
codes. This Grobner bases setting turns out to be strongly related to the 
combinatorics of the codes. 
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■ Introduction 

00 

. It is well known that the complexity of Grobner bases computation heavily 

' depends on the term orderings, moreover, elimination orderings often yield a 

. greater complexity. This remark led to the so called FGLM convertion pro- 

' blem, i.e., given a Grobner basis with respect to a certain term ordering, find 

, a Grobner basis of the same ideal with respect to another term ordering. One of 

the efficient approaches for solving this problem, in the zero-dimensional case, 
is the FGLM algorithm (see [TT]L 

The key ideas of this algorithm were successfully generalized in ^21 with the 
objective of computing Grobner bases of zero-dimensional ideals that are de- 
termined by functionals. In fact, the pioneer work of FGLM and was the 
Buchberger-Moller's paper (cf. Authors of P used the approach of [5j and 
some ideas of ^T] for an efficient algorithm to zero-dimensional schemes in both 
affme and projective spaces. In 4 similar ideas of using a generalized FGLM 
algorithm as a pattern algorithm were presented in order to compute Grobner 
basis of ideals of free finitely generated algebras. In particular, it is introduced 
the pattern algorithm for monoid and group algebras . In [3] a more general 
pattern algorithm which works on modules is introduced, many things behind 
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of this idea of using linear algebra are formalized, notions like "Grobner tech- 
nology" and "Grobner representations" are used. There are other approches 
which also generalized similar ideas to some settings, behind of all these works 
is the essential fact of using linear algebra techniques to compute in "Grobner 
bases schemes" . 

The main goal of this paper is to show the application of techniques to 
linear codes like the ones in FGLM and subsequent works, which comes from 
an specification of the pattern algorithm for monoid algebras given in 0], i.e. 
by taking an algebra associated to a linear code. 



1 Preliminaries 

The case of the algebra associated to a linear code we are going to introduce 
is connected with an ideal of a free commutative algebra; therefore, we will 
restric ourselves to the formulation of a pattern algorithm for a free commutative 
algebra. Similar settings can be performed in a free associated algebra or over 
modules (see 

Let X :— {xi, . . . , x n } be a finite set of variables, [X] the free commutative 
monoid on X, K a field, I an ideal of K[X], 1(F) the ideal of K[X] generated 
by F C lsr[-X"], K[X\/I the residue class algebra of K[X] module I. Let us 
denote by 1 the empty word in [X], L(u) the length of the word u 6 [X], and 
Gard(C) the cardinal of the set C. Let now -< be a semigroup total well ordering 
on [X] (such an ordering is also called admissible), then for / G K[X] \ {0}, 
T^(/) is the maximal term of / with respect to -<, LC^(f) is the leading 
coefficient of / with respect to -<. Similarly, for F C A"[X], T^{F} is the set 
of maximal terms of non-zero polynomials in F, T^(F) is the semigroup ideal 
generated by T{F}. Moreover, for the sake of simplicity in notation, U^(F) 
will be used instead of U(T^(F)), where U lies in {G, N, B, I}. Of course, given 
an ideal I and two different admissible orderings -<i and -<2, in general we have 
U(T^ 1 (Tj) 7^ U(T^. 2 (I)). Notwithstanding this strong dependency on -<, while 
a single admissible ordering -< is considered, so that no confusion arise, we will 
often simply write U(F) for U^(F). 

Let r C [X] be a semigroup ideal of [X], i.e., for u G [X] and t G r, tu G r. 
Then, it is well known that r has a unique subset G(t) of irredundant generators 
(probably infinite). In the case of / a zero-dimensional ideal, for r = T(I), G(t) 
is always finite. We are going to introduce for r some notation and terminology, 
which are similar to those introduced in |12| . 

Pred(w) := {u G [X] \ 3 x G X (w = ux)} (the set of predecessors of w), 
N(t) := {s G [X] s £ t} (outside o/r), 

B(t) :— {w G r Pred(w) C N(t)} (border o/r), 

7(r) := t \ B(t) (interior o/r). 

We remark that w G r lies in G(r) if and only if all its proper divisors are 
in N(t) (that is if Pred(w) C N(t)). In the following proposition, some basic 
results concerning r and its regions are summarized. Although they are very 
easy to prove, their importance is crucial for FGLM techniques. 
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Proposition 1 (Properties of the semigroup ideal regions). i. For each 
w E t there exist u G [X] and v € B(r) s.t. w = vu. 

ii. For x G X : 

(a) Ifue N(t), then ux G N(t) U B(t). 

(b) Ifue B(t), then ux G B(t) U I(t). 

(c) If u G I(t), then ux G I{t). 

Hi. N(t),N(t) U G(t), N(t) U B(t) are order ideals, i.e., if u belongs to one 
of these subsets and v divides u, then v also belongs to the corresponding 
sets. 

Theorem 1 (The vector space of canonical forms modulo an ideal). Let 

Spanx(N ^(1)) be the K -vector space whose basis is 7V_<(J). Then the following 
holds: 

i. K(X) = I © SpanK{N^{I)) (this sum is considered as a direct sum of 
vector spaces). 

ii. For each f G ■K'l-X"] there is a unique polynomial of SpariK{N^(I)), de- 
noted by Can(f, I, -<) such that f — Can(f, I, -<) G /; moreover: 

( a) Can(f, /,-<) = Can(g, I, -<) if and only if f — g G I. 

(b) Can(f,I,~<) — if and only if f G I. 

Hi. There is a K -vector space isomorphism between K[X]/I and Spanx(N^ (I)) 
(the isomorphism associates the class of f modulo I with the canonical 
form Can(f, I, -<) of f modulo I). 

Can(f,I,~<) is called the canonical form of / modulo /. We use simply 
Can(f,I) if the ordering used is clear from the context. 

We assume the readers to be familiar with definition and properties of 
Grobner bases (see [2] for an easy to read introduction to Grobner bases). 

Proposition 2 (Characterization of zero-dimensional ideals). Let G be 

a Grobner basis of I with respect to -i. Then, I is a zero- dimensional ideal (i.e. 
dirriKK[X]/ 1 < oo) if and only if N^(G) is finite. Moreover, in such a case, 
dim K K[X]/I = Card(N^(G)). 

Definition 1 (Border basis). The border basis of I with respect to -< is the 
subset B(I, -<) G I defined by: 

B(I, -<) := {w - Can(w, I,-<) | w G B^(I)} (the B-basis of I). 

Note that the £>-basis of J is a Grobner basis of / that contains the reduced 
Grobner basis. 
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1.1 Matphi matrices and Grobner representation 

The word Matphi appears by the first time in JT] to denote a procedure that 
computes a set of matrices (called matphi matrices) s.t. there is one matrix for 
each variable in X and they describe the multiplication structure of the quotient 
algebra K[X]/I, where 7 is a zero dimensional ideal. We often refer to this set 
of matrices as the matphi structure. 

Definition 2 (Grobner representation, Matphi structure). Let I be a 

zero-dimensional ideal of K[X], lets = dim(K[X]/I). A Grobner representation 
of I is a pair (N, <f) consisting of 

i. N — {N x , ...,N S } s.t. K[X]/I = Span K (N), and 

ii. cf> :— {(f>(k) | 1 < k < n}, where (^(k) are the square matrices (f>(k) :— 
(o^-)ij s - 1 - f or alll <i<s, NiX k ee/ a'-jNj. 

<f> is called the matphi structure and the <j){k) 's the matphi matrices. 

See |2] for a more general treatment of these concepts. Note that the matphi 
structure is indepent of the particular set N of representative elements of the 
quotient K[X]/I. In addition, the matphi matrices allow to obtain the class of 
any product of the form NiXk as a combination of the representative elements 
(i.e. as a linear combination of the basis N for the vector space K[X]/I). 



In this section we present a generalization of the FGLM algorithm for free 
commutative algebras, which allows to solve many different problems and not 
only the clasic FGLM convertion problem. The procedure we are presenting 
is based on a sort of black box pattern: in fact, the description of the steps 
5 and 6 is only made in terms of their input and output. More precisely, we 
are assuming that a term ordering -<i is fixed on [X], I is a zero-dimensional 
ideal (without this restriction the algorithm does not terminate), and that the 
K- vector space SpanK(N^ 1 (I)) is represented by giving 

• a TT-vector space E which is endowed of an effective function 



which, for each finite set {v\, . . . , v r } C E of linearly independent vectors 
and for each vector v € E, returns the value defined by 



2 The FGLM pattern algorithm 



LinearDependency[w, {vi, . . . , v r } 




if V = YJi=l 

if v is not a linear combination of {v\, . . . , v r }. 



• an injective morphism £ : Spanx{N ^ t (I)) i— > E. 
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This informal approach allows a free choice of a suitable representation of the 
space SpariK(N -<i(-0) regarding an efficient implementation of these techniques 
and a better complexity. Moreover, as an aside effect, it enables us to present 
this generalization in such a way that it can be applied on several more particular 
patterns and helps to make key ideas behind the FGLM algorithm easier to 
understand. Let us start making some references to some subroutines of the 
algorithm. 

Insert Next s[w, List, -<] inserts properly the products wx (for x £ X) in List 
and sorts it by increasing ordering with respect to the ordering The reader 
should remark that InsertNexts could count the number of times that an 
element w is inserted in List, so w £ N^(I) UT^{G} if and only if this number 
coincide with the number of variables in the support of w, otherwise, it means 
that w £ T^(I) \ T^{G}, see [TT], this criteria can be used to know the boolean 
value of the test condition in Step 4 of the Algorithm 
NextTerm[List] removes the first element from List and returns it. 

Algorithm 1 (FGLM pattern algorithm). 

Input: -<2, o, term ordering on [X]; £ : Spanx(N^ <1 (I)) i— ► E. 

Output: rGb(I,~^2), the reduced Grobner basis of I w.r.t. the ordering -<2- 

1. G := 0; List := {1}; N := 0; r := 0; 

2. While List ^ do 

3. w := NextTerm[List]; 

4. If w T^ 2 (G) (if w is not a multiple of any element in G) then 

5. v := ^(Can(w,I,<i)); 

6. A := LinearDependency[u, {v\, . . . , v r }}; 

7. IfA^ False then G := GU{w-J2l =1 (where A = (Ai, . . . , X r )) 

8. else r := r + 1; 

9. v r :— v; 

10. w r := w; N := N U {w r }; 

11. List := InsertNexts [uv, List, -<%]; 

12. Return[G]. 

Remark 1. i. A key idea in algorithms like FGLM is to use the relationship 
between membership to an ideal I and linear dependency modulo I, namely 
Vc t £ K, Si e K[X}: 

Y^i=i c i s i ^ I \ {0} {' s i> • • ■ i s r} * s linearly dependent modulo I. 
This connection with linear algebra was used for the firts time in Grobner 
bases theory since the very begining (see JBl). 

ii. Since each element of N^ 2 (I) U B^ 2 (I) belongs to List at some moments 
of the algorithm and List C N^ 2 (I) U B^ 2 (I) at each iteration of the 
algorithm, it is clear that one can compute B{I , -<2) or the Grobner re- 
presentation (N^ 2 (I), <fi) of I just by eliminating Step 4 of the algorithm 
and doing from Step 5 to Step 11 with very little changes in order to built 
those structures instead of rGb{I , -<2) • 
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Hi. Note that Step 5 and 6 depends on the particular setting. In Step 5 it is 
necessary to have a way of computing Can(w, I, -<i) and the corresponding 
element in E, while in Step 6 we need an effective method to decide linear 
dependency. 

iv. Complexity analysis of this pattern algorithm can be found in for the 
more general case of free associative algebras, and for a more general 
setting in ' L 3L Of course, having a pattern algorithm as a model, it is 
expected that for particular applications, one could do modification and 
specification of the steps in order to improve the speed and decrease the 
complexity of the algorithm by taking advantage of the particular structures 
involved. 

2.1 The change of orderings: a particular case 

Suppose we have an initial ordering -<i and the reduced Grobner basis of / for 
this ordering, now we want to compute by the FGLM algorithm the new reduced 
Grobner basis for a new ordering -<2- Then the vector space E is K s , where 
s = dim(K[X]/ 1). In Step 5, Can(w, I, -<i) can be computed using the reduced 
Grobner basis rGb(I, -<i) and the coefficients of this canonical form build the 
vector of E corresponding to this element (the image by the morphism £). Then 
Step 6 is perfomed using pure linear algebra. 

3 FGLM algorithm for monoid rings 

The pattern algorithm is presented in 4 for the free monoid algebra, we will 
restrict here to the commutative case. Let M be a finite commutative monoid 
generated by g\, . . . , g n ; £ : [X] — > M, the canonical morphism that sends Xi to 
g%; a C [X] x [X], a presentation of M defined by £ (a = {(w, v) \ £(w) = £(v)}). 
Then, it is known that the monoid ring JTfM] is isomorphic to K[X]/I(a), 
where 1(a) is the ideal generated by P(a) — {w — v | (w,v) £ cr}; moreover, 
any Grobner basis G of 1(a) is also formed by binomials of the above form. In 
addition, it can be proved that {(w,v) \ w — v € G} is another presentation of 
M. 

Note that M is finite if and only if I = 1(a) is zero-dimensional. We will show 
that in order to compute rGb(I), the border basis or the Grobner representation 
of /, one only needs to have M given by a concrete representation that allows 
the user to multiply words on its generators; for instance: M may be given by 
permutations, matrices over a finite field, or by a more abstract way (a complete 
or convergent presentation). Accordingly, we are going to do the necessary 
modifications on Algorithm \l\ior this case. 

We should remark that in this case -<i=^2, then at the begining of the 
algorithm the set N^ ± (I) is unkown (which is not the case of the change of 
orderings). It could be precisely a goal of the algorithm to compute a set of 
representative elements for the quotient algebra. 
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Now consider the natural extension of £ to an algebra morphism (£ : K[X] i— > 
if[M]), note that the restriction of £ to Span K (N^ 1 (/))(£: Span K {N^ 1 (I)) i-> 
if[M]) is an injective morphism; moreover, = £(Cem(iy, 7, -<i)), for all 

w G [X], Therefore, the image of Can(w,I, -<i) can be computed as £(w), and 
the linear dependency checking will find out whether w is a new canonical form 
(i.e. w 6 AT Xl (i")) or not (i.e. to £ T Xl (rG6(J, ^1))). Hence, Step 5 will be 
v := £(u)gi, where it £ Pred(w) and mi, = 

Moreover, let u>i, . . . , w r be elements of N^ il (i) and Vi = £(it>i), for 1 < i < r. 
Then LinearDependency[v, {v\, . . . , v r }) can be computed as 

{vj if v = Vj, for some j £ 

False otherwise. 
Finally, Step 7 changes into: 

If A ^ False then G := G U {w - Wj }. 

Remark 2. i. This example shows that the capability of the K -vector space 
E w.r.t. Linear-Dependency, that is demanded in the Algorithm 1, is 
required only on those sets of vectors {v\, . . . , v r , v} that are built in the 
algorithm, which means in this case that LinearDependency is reduced 
to the Member checking, i.e., v is linear dependent of {vi, . . . ,v r } if and 
only if it belongs to this set. 

ii. When a word w is analyzed by the algorithm, all the elements in Pred(w) 
have been already analyzed (£(u) is known for any u G Pred(w) ), this is the 
case whenever -<i is an admissible ordering. Therefore, the computation 
of £(u>) is immediate. 

We will show the case of linear codes as a concrete setting for an application 
of the FGLM pattern algorithm for monoid rings, where the monoid is given by 
a set of generators and a way of multiply them. 

4 FGLM algorithm for linear codes 

For the sake of simplicity we will stay in the case of binary linear codes, where 
more powerfull structures for applications are obtainned as an output of the 
corresponding FGLM algorithm (for a general setting see From now on 

we will refer to linear codes simply as codes. 

Let F2 be the finite field with 2 elements. Let C be a binary code of dimension 
k and length n (k < n), so that the n x (n — k) matrix H is a parity check matrix 
(c • H = if and only if c G C). Let d be the minimum distance of the code, 
and t the error-correcting capability of the code (t = [^j^l , where [x] denotes 
the greater integer less than x). Let B(C,t) = {y G F£ | 3 c G C (d(c,y) < f)}, 
it is well known that the equation eH = yH has a unique solution e with 
weight(e) < t, for y G B(C,t). 

Let us consider the free commutative monoid [X] generated by the n va- 
riables X := {xi, . . . ,x n }- We have the following map from X to ip '■ 
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X — > Fj, where Xj i— > (the z-th coordinate vector). The map ip can be ex- 
tended in a natural way to a morphism from [X] onto F£ , where VKIliLi x i) = 
(/3imod2,...,/3 n mod2). 

A binary code C defines an equivalence relation Rq in given by (x, y) £ Rc 
if and only if x — y e C. If we define £(u) := ip(u)H, where u e [X], the above 
congruence can be translated to [X] by the morphism ip as u =c w if and only 
if (ip(u),ip(w)) € i?c, that is, if = The morphism ^ represents the 

transition of the syndromes from F£ to [X]; therefore, £(u>) is the "syndrome" 
of w, which is equal to the syndrome of ip{w). 

Definition 3 (The ideal associated with a binary code). Let C be a binary 
code. The ideal 1(C) associated with C is 

1(C) := ({w-u\ eH-eH})civ[X]. 

5 The algorithm for binary codes 

The monoid M is set to be ¥2~ k (where the syndromes belong to). Doing 
9i ■= £(xi), note that M = F2~ fe = (g±, . . . ,g n ). Moreover, a := Rc, hence 
1(a) = 1(C). Let -< be an admissible ordering. Then the FGLM algorithm 
for linear codes can be used to compute the reduced Grobner basis, the border 
basis, or the Grobner representation for -<. 

Algorithm 2 (FGLM for binary codes). 

Input: n,H the parameters for a given binary code, -< an admissible ordering. 
Output: rGb(I(C),~<). 



1. 


List := {1}, N := 0, r := 0, G = {}; 


2. 


While List ^ do 


3. 


w := NextTerm [List]; 


4. 


Ifw£T(G); 


5. 


v := £(w); 


6. 


A := Member [v, {vi, . . . , v r }]; 


7. 


If A ^ False then G := G U {w - wj}; 


8. 


else r := r + 1; 


9. 


v r := v' ; 


10. 


w r := w, N := N U {w r }; 


11. 


List := InsertNext[u> r , List}; 



12. Return [G]. 
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In many cases of FGLM applications a good choice of the ordering -< is a 
crucial point in order to solve a particular problem. In the following theorem 
it is shown the importance of using a total degree compatible ordering (for 
example the Degree Reverse Lexicographic). Let us denote by <t a total degree 
compatible ordering. 

Theorem 2 (Canonical forms of the vectors in B(C,t)). Let C be a code 
and let Gt be the reduced Grober basis with respect to <t- If w £ [X] satisfies 
weight(-0(Can(w, Gt))) < t then ?p(Can(w,GT)) is the error vector correspon- 
ding to ip(w). On the other hand, if weight (il>(Can(w, Gt))) > t then ip(w) 
contains more than t errors. 

Proof. If we assume that weight ( "0( Can (w, Gt))) < t then, we can infer at 
once that ip(w) <E B(C,t) and ip(Can(w,GT)) is its error vector (notice that 
£{ w ) = £(Can(w, Gt)) and the unicity of the error vector). 

Now, if weighty (Can (u>, Gt))) > t, we have to prove that ifj(w) ^ B(C,t). 
It is equivalent to show that weight(^(Can(ui, Gt))) < t if %p(w) £ B(C,t). Let 
tp(w) be an element of B(C,t) and let e be its error vector then, weight(e) < t. 
Let w e be the squarefree representation of e. Note that weight(e) coincides 
with the total degree of w e ; accordingly, L(w e ) < t. On the other hand, 
Can(w,GT) <t w e , which implies that L(Can(w,GT)) < £(w e ) (because <t 
is degree compatible). Hence, weight (ip(C an (w,Gt))) < L(Can(w,GT)) < t. 
□ □ 

The computation of the error-correcting cability of the code t can be done 
in the computing process of Algorithm [21 (see the example in Section 15.11 and 
|7| )■ The previous theorem allows us to use the computed reduced Grobncr 
basis for solving the decoding problem in general binary codes, but also with 
such a powerful tool available, it is expected to be able to study the structure 
of the codes, like some combinatorics properties. Some possible examples are 
the permutation-equivalence of codes (see [5]), and some problems related with 
binary codes associated with the set of cycles in a graph (finding the set of 
minimal cycles and a minimal cycle basis of the cycles of a graph), see 

To generalize Theorem |2 for non binary linear codes have some conflicts 
with the needed ordering; however, the FGLM algorithm can be still used to 
compute the border basis or a Grobner representation for the ideal 1(C) and 
it will be possible to solve the problems that one can solve with the reduced 
Grobner basis in the case of binary codes. Those problems are explained in [3- 
In addition, 5 contains some results and examples about the application of this 
setting to general linear codes and, in binary codes, for studying the problems 
of decoding and the permutation-equivalence. 

5.1 An example 

Let C be the linear code over Fj] determined by the parity check matrix H given 
below. The set C of codewords is given in the right hand side. The minimum 
distance is d = 3, so, t = 1, the numbers of variables is 6, <t is set to be the 
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Degree Reverse Lexicographic ordering with Xi+i >t %i 
of the computation will be described. 



Only essential parts 



H = 



C = {(0, 0, 0, 0, 0, 0), (1, 0, 1, 1, 0, 0), 
(1, 1, 1, 0, 0, 1), (1, 1, 0, 0, 1, 0),(0, 1, 0, 1, 0, 1), 
(0, 1, 1, 1, 1,0), (0, 0, 1,0, 1, 1), (1, 0, 0, 1, 1, 1)}. 



Application of Algorithm |2J List := {1}; N := {}; r := 0; w := 1; £(1) = 
(0,0,0); N := AU{1} = {1}; £(iV) := {(0,0,0)}; List := { Xl ,x 2 ,x 3 ,x 4 ,x 5 ,x 6 }; 
w:=x i; £(asi) = (1, 1, 1); N := {1,^}; £(JV) := {(0, 0,0), (1, 1, 1)}; 
After analyzing xq we are at the following stage: 

N := {1,xi,X2,X3,X4,X5,xq}, and List — {x\, X\x 2 , XiXs, X1X4,, Xxx$, X\X6,x 2 , 
X2X3, x 2 X4, x 2 x 5 , x 2 x 6 , x% x 3 x 5 , x 3 x 6 , x%, XiX bl X4X 6 , x\, x 5 x 6 , x%}. 

There is still one element left in N because there are 7 elements in N of 
a total of 8 (2 6 ~ 3 ). Taking the elements of List from x\ to X1X5 they are a 
linear combination of elements already in N (their syndromes are in the list of 
syndromes computed £ (N)). Therefore, G :— {xf — 1, X\x 2 — x§, X1X3— 14, X1X4— 
%3) X1X5 — x 2 }, for example x\x 2 — x$ is obtained, bacause when w — x\x 2 , first 
note that Pred(w) C N, which means that it is either a new irreducible element 
or a head of a binomial of the reduced basis. Then £(x\X 2 ) is computed and 
we got that (,{x\X 2 ) = £(#5) . This means that x\x 2 — x$ belongs to G. Also 
x±x 2 is the first minimal representation which is not in N, this implies that 
t = weight (V* (2; 1 £2)) — 1 (see 0). The next element in List, w = x±Xq, is the 
last element that will be included in N and the corresponding multiples will be 
included in List. 

From this point, the algorithm will just take elements from List and it ana- 
lyzes in each case whether it is in T{rGb(I(C), <t)} (like X2X3) or in 
T(rGb(I(C), <t) \ T{rGb(I(C), <t)} (like xix 2 Xe), this process is executed un- 
til the List is empty when the last element xix^ of the list is analyzed. Finally 
the reduced Grobner basis for <t is 

G := {x\ - \,xxx 2 — £5, £123 - x 4 ,xia;4 — x 3 ,xix 5 - x 2 ,x\ - ^,X 2 x 3 — 
■Xx,X 2 Xq — X4,x1-1, £3X4 - Xi , £3X5 - x e , x 3 x 6 - x 5 , x\ - 



XiXq, X 2 X4—X 6 , X 2 X 5 

1, £4X5 - xix e , x 6 X4 - x 2 , xl - 1, x 5 x 6 - x s , x% — 1}. 

Now let us assume that a vector y = (1, 1, 1, 0, 1, 0) is received, the corres- 
ponding word is w = x\x 2 x 3 X5. Then we compute w e = Can(w,G) = 23. As 
weight (ip(w e )) = 1 (t = 1 is the error-correcting capability); therefore, the error 
vector is e = (0, 0, 1, 0, 0, 0), and the codeword is c = (1, 1, 0, 0, 1, 0). 
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